Contribution aux Journees Relativistes 1979, ed. I. Moret-Bailly & C. Latremoliere, 
Anneau d'accretion sur les trous noirs" pp 166-182 (Faculte des Sciences, Angers, 1979). 



GENERALISED MASS VARIATION FORMULA 
FOR A STATIONARY AXISYMMETRIC STAR OR 
BLACK HOLE WITH SURROUNDING ACCRETION DISC 

Brandon Carter 

(Nj \ G.A.R., Observatoire de Paris, 92 Meudon. 

O ■ 
O ' 

<n : 

Abstract The general relativistic mass-energy variation formula for axisymmetric equi- 
iJ^ \ librium states of a self gravitating system is developed in the particular case for which the 

^ ■ relevant matter consists of a perfectly conducting multiconstituent fluid (or superfluid) 

with stationary circular or convective motion. 

> '■ 

co : 

t— i ■ 1. Introduction 

^° '■ 

This communication presents an extension of the study 1 of mass variation for- 
mulae for stationary axisymmetric systems - relativistic stars and black holes with 
surrounding material discs - that is contained in my contribution (chapter 6) to the 
Einstein Centenary Survey: General Relativity, ed S.W. Hawking and W. Israel (Cam- 
bridge U.P. 1979); all numbered references here refer to equations presented therein. 
£l ', The present work starts from the general formula 

a ■ 

8M=^-8 I G a b k a d£ b - / G cd h cd k a d£ a + fi H 5J U + ^5A (1) 

■ 57T J 07T J 07T 

rS ' 

(equation (6.342) 1 ) for the variation of the total asymptotically defined mass M of the 
system between two neighbouring stationary axisymmetric equilibrium states, where 
are the angular velocity, angular momentum, decay parameter ("surface 
gravity") and surface area of the central black hole - if it is present - while G a b is 
the Einstein tensor, h c d the variation of the metric tensor, k a the stationary Killing 
vector, and the integration is taken over a spacelike hypersurface £ extending out to 
the asymptotically flat region with inner boundary on the black hole if it is present. 

2. Electromagnetic and material contributions 

We now introduce the energy momentum tensor T a b and the corresponding mo- 
mentum vector II a = —T a b k b (decomposed into electromagnetic and purely material 
parts, labelled by suffices p and m ) as follows. When we substitute from the Einstein 
equations into the fundamental kinematic mass variation formula (1) we obtain 

5M - fi H 5J H - 5 A = 5 J T a b k b d£ a - \ J T b c h c b k a d£ a (2) 
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i.e. 

6M-n*6J s --^6A+6j n F a {k} d£ a + \ J T b c h c b k a dE a 

= ~ S J n M a i k } f 1 ^" dE « ( 3 ) 

To evaluate the electromagnetic contribution we first use the analogue of (6.306) 
for the momentum flux of a stationary field, namely 

-77/ = A F k a + (k b A b )j a + V b (A c k c D F ba ) (4) 

and then work out 

\g\- 1 / 2 5{\g\ 1 / 2 K ¥ k a ) = \(T r b e h\- D bc 8F bc )k a 

= {\T ¥ bc h hc - j b 5A b )k a + 2V b {D/ b 5A c k^) , (5) 

where the electromagnetic action density and displacement tensor are given - by (6.193) 
and (6.194) 1 - as A F = (l/167r)F a6 F ba and D F ab = (l/4yr)F ab . 

With the aid of the Green theorem (and using the asymptotic flatness conditions 
to eliminate a surface contribution at infinity) we are thus able to obtain 

5 I n F a d£ a -i f T F b c h c b k a dX a = j k b A b 5{fdX a ) + ± I fc[yi(M 6 )<ffi 

-\ £ \k c A c D ba dS ab - £ D/ b (8A c )k a l dS ab (6) 



3. Perfectly conducting fluid 

The result so far does not depend on the particular field equations satisfied by any 
matter that is present. In order to show the relationship with the first law of ordinary 
thermodynamics we shall now consider the explicit form of the right hand side of (3) 
in a particularly ideal case namely that of a perfectly conducting multiconstituent fluid 
medium. Previous descriptions of the first law by Bardeen 1973 2 , Bardeen, Carter and 
Hawking 1973 3 have been restricted to the more specialised case of a non-conducting 
medium, electromagnetic effects being taken into account in the more extended treat- 
ment given by Carter 1973b 4 . 

The more general treatment presented here is based on the relativistic theory of 
a multiconstituent perfect fluid obtained as the natural extension of the two constituent 
case discussed by Carter 1976 5 . According to this theory (which includes the Landau- 
London two-fluid model for superfluidity as a special case in the non-relativistic limit) 
the material equations of motion are derivable from a Lagrangian scalar A M which is a 
covariant function of a set of conserved current vectors n x representing the fluxes of a 
set of independent constituents labelled by an index x, one of the constituents being 
usually (except in zero-temperature models) the entropy. (Thus in a typical application 
x might take the values i, i, n with n 1 and n a representing fluxes of positively charged 
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ions and negatively charged electrons respectively, while n m is identified with an entropy 
flux vector s. 

Under the influence of arbitrary infinitesimal variations of the fluxes n and of the 
metric g the function A M will have an infinitesimal variation of the form 

5K=^Jn a x + \^ a n b x h%. (7) 

The covectors \x x so defined may appropriately be described as the chemical 'potential 
one-forms of the corresponding conjugate constituents. (In particular the covector con- 
jugate to the entropy flux vector s may appropriately be denoted by © and interpreted 
as a temperature form.) Their constituent indices are written upstairs in view of the 
fact that they transform contravariantly with respect to the fluxes under a change of 
chemical basis in the vector space of constituents. 

When the medium interacts with an electromagnetic field the total Lagrangian 
will be taken to have the simple additive form 

A = A M + A a j a + A F (8) 

where A F is the electromagnetic Lagrangian (given by (6.193) 1 ) and where the electro- 
magnetic current j a appearing in the interaction term is given by 

j = e x n x (9) 

where e x represents the mean electric charge per particle of the xth constituent. Like 
the chemical potentials, the charges behave as components of a contravariant vector in 
the chemical constituent space. (It would in fact be easy to work with a more general 
Lagrangian involving a more intimate relationship between the fluxes n x and the field F_. 
The restriction to the simple additive form (8) - which means that polarisability effects 
are excluded - has been made merely to facilitate the interpretation of the formalism.) 
Variation of the total Lagrangian now gives 

SA = iv x dn« + Cf - j«)5A a - V c (ZT 6 <L4 b ) 

+ J(n>£ - D™V b A c + V c {D™A b ) + (tf - j a )A h )h a b (10) 
where we have used the abbreviation 

j; = v b Df (ii) 

in accordance with the scheme introduced in section 6.4 1 , and where we introduce the 
notation 

n x = + e x A (12) 

for what will play the role of an effective four-momentum per particle for the correspond- 
ing constituent. If we required the volume integral of A to be invariant under arbitrary 
variations of the fields A and n x we would obtain not only the correct equations of 
motion 

l=J (13) 
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for the electromagnetic field, but also the too restrictive requirement that the momenta 
7r x vanish altogether. To obtain reasonable equations of motion for the chemical con- 
stituents one must demand invariance not under arbitrary current variations but only 
under variations of the restricted class that arise from infinitesimal displacements of 
the world lines of the idealised particles, or directly from variations of g a \>. Now it can 
easily be seen (e.g using the formulae given by Carter 1973a 6 ) that an infinitesimal 
displacement £ x of the world lines of a conserved flux of particles (acting in conjunction 
with a variation h of the metric) induces a corresponding variation 

Sn a = -\n a h c c - n a V b ^ b + n b V b C - ^ b V b n a (14) 

in the corresponding current vector n. 

The equations (7) to (13) were chemically covariant (i.e. unaffected by changes 
in the chemical basis of the space of constituents) in form although not necessarily in 
content (since the specific representation of the Lagrangian will in general depend on the 
chemical basis). However the choice of a particular set of chemical constituents implies 
the choice of a preferred basis in the vector space of constituents, so that the system of 
equations that results will not have even the appearance of chemical covariance. It is also 
to be remarked that although the identities derived in section 6.4. 1 1 (which depended 
only on the covariance of the Lagrangian under general co-ordinate transformations of 
the space-time manifold as opposed to the chemical constituent vector space) will still be 
formally valid for the material lagrangian A M - indeed (6.135) 1 has been implicitely used 
in the derivation of the form (7) - and for the total Lagrangian A, the quantities involved 
cannot be given the same direct interpretation (as representing energy-momentum etc. 
...) for a restrained variational theory as in the free field theory for which they were 
originally intended. In the present theory the appropriate energy-momentum tensor 
cannot be obtained directly from the analogue of (6.138) 1 but can be read out from the 
expression 

M- 1/2 5(M 1/2 A) = T ab h ab + (f-j-)SA a - + V a (2nfn^S ~ A b D* b ) (15) 

where the £ x (which do not behave as components of a well behaved vector under 
chemical basis transformations) are a set of arbitrary displacement vector fields acting 
on the world lines of the corresponding fluxes, and where the fields f x (to which the 
same remark applies) are to be interpreted as generalised 4-force densities acting on the 
particles of the corresponding constituents. One finds that the force covectors have the 
explicit form 

l x =n }xl -dK x +K x V-n lx{ (16) 

where the inclusion of the chemical indices within straight bars is introduced to denote 
suspension of the summation convention. 

Each of the variational field equations 

r = o (i6) 

can be decomposed by first contracting with the corresponding flux n x , to give first the 
conservation law 

V • n x = 18 
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for each constituent, and then the canonical equation of motion 



n m ■ drr* = (19) 

which determines the acceleration of the corresponding world lines. This last equation 
may be decomposed into a material and an electromagnetic part in the form 

n lxl -dn x =F-j x (20) 

where 

J x = e x n }x} (21) 

is the electric current density associated with the xth constituent. It is to be remarked 
that the condition (19) is precisely what is required for the flux of generalised vorticity 

w x = dn x (22) 

to be conserved along the flow lines in the strong sense discussed in section 6.4.5 1 , 
meaning that 

(an lx )£w x = (23) 

for an arbitrary scalar field a, which implies that the flux w x over any 2-surface, which 
by Stokes theorem is equal to the circulation 



C x = j) K x ■ dx (24) 

around its boundary circuit c say, is conserved in the sense of being unchanged when 
the circuit (or the original 2-surface enclosed by it) is dragged at an arbitrary rate along 
the flow lines. 

Unlike the field equations themselves, the expression for the energy momentum 
that is read out from (15) has a chemically covariant form. After subtracting out the 
Maxwellian contribution in accordance with the decomposition (6.299) 1 we are left with 
a material contribution to the energy-momentum given by 

T Mb = ™>l ~ *»9% (25) 

where \& M is a thermodynamic potential density that is related to the material Lagrangian 
by 

* M =n>--A M . (26) 

Provided that the electromagnetic field equations (13) are satisfied, the expression 
(6.299) 1 of the total metric energy momentum tensor as the sum of metric contributions 
T F ab and T M ab can be replaced by an expression in terms of canonical contributions in 
the form 

T a b = T M a cb + T\ + V c {D™A b ) (27) 

where T a b are components of the canonical electromagnetic energy momentum tensor 
introduced in (6.145) 1 , and where the corresponding canonial material energy momen- 
tum tensor is defined by 

T M^ = «X-V fc - (28) 
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Subject to the field equations, the total energy momentum tensor will of course be 
preserved, i.e. 

V 6 T b a = (29) 
and hence the separate parts will satisfy 

V b T M b a = F ab j b = -V b T F b a , (30) 
V b T M b ca =j b V a A b = -V b T b a . (31) 



For any vector field £ we define associated canonical and metric energy - momen- 
tum flux vectors analogous to those introduced in section 6.4 1 by 



n{0 = n M {0 + n F {0 , n M a {0 = -T m \e , (33) 

the total being related by 

V a {i} = n a {i} + V b (D ba A c C) . (34) 
It is obvious from (29) that the total metric energy - momentum fluxes will be conserved 

— * 

whenever £ is a Killing vector, and hence that the same applies to the canonical energy 
- momentum flux, i.e. 



££<7 = 



v • n{£} = (35) 
V • V{£} = . (36) 



However the canonical energy - momentum flux also has the property that is material 
part will be conserved separately whenever the fields are themselves invariant. Clearly 
the equation of motion (19) implies (by the Cartan formula (6.20) 1 ) that we shall have 
a generalised Bernoulli theorem to the effect that 

££tt x =0 => n lxl -V(f 1L x ) = , 

(37) 

i.e. that £ • 7r x is a constant of the motion, and hence, using (18) we see that 



6£tl- = o , ^ M = o 

&g = 



V-P MC {£} = 0. (38) 



This is to be compared with the corresponding result (6.157) 1 for the electromagnetic 
component, which tells us that the same result would follow from purely electromagnetic 
invariance, i.e. 

(£A = 



&9 = 



V-P MC {£} = 0. (39) 



It may be noticed that across any invariant hypersurface, i.e. one whose surface element 
satisfies 

£ ■ dE = (40) 
6 



the canonical material energy - momentum flux is interpretable as the total rate of 
transport of the Bernoulli quantity, i.e. 

P M « c {£}dE a = (£- K x)dN x (41) 

where dN x denotes the number flux of the xth constituent across dE , i.e. 

dN x = -n x ■ dE . (42) 



4. Global energy and angular momentum 

If we now consider the particular case of a stationary and/or axisymmetric system, 
in which there will be a conserved energy and/or angular momentum per particle given 
by 

S x = -k ■ n x (43) 
J* = m ■ n x (44) 

then we shall have 

Kc^^a = ~^N X = "£Z c l ■ (45) 



X 



C M dS « = J x dN x =-J2 J C M ■ d ^ ( 46 ) 



X 



across any respectively stationary or axisymmetric hypersurface, and in particular across 
a black hole horizon, where we define the fluxes 

4 X c=^,*, (47) 

J M x c =J x n ]x] (48) 
which evidently satisfy the separate conservation laws 

V ' 4 X c = (49) 

V'^=0. (50) 

In the case of angular momentum, the total exterior source contribution given by 
(6.322) 1 may be expressed by 



E / 4 • (51) 

x J * 



for any axisymmetric hypersurface E, but the fact that E could not be invariant under 
k prevents us from obtaining an equaly simple formula for M MC . 

We now proceed from (6) as in Section 6.6.1 using the strict stationary boundary 
conditions (6.314) 1 to evaluate the surface integrals on the horizon in the form 

5 I \k c A c D b F a dS- <f Df(8A c )k a ^ dS ab = $ H 5Q H + H 5J FH (52) 
Js Js 



We next perform the analogous steps in the evaluation of the material contributions, 
using the formulae of the previous section which enable us to write firstly 



~n M a {k} = (A M - n b x y,% )k a + (fcVfK (53) 

and 

|$r 1/2 *(|0| 1/2 A M ) = ±(T M b c hl + n b x » x b h c c ) + ^ x 5n b . (54) 

Since we now only have source terms to deal with we proceed directly (i.e. without 
recourse to Green's theorem) using (6.399) 1 to obtain 

-5 I n^ a -\ ! T m b c hlk a dZ a = ! fcVS*(n-dE a ) + i f fc[» n -](^)dE (55) 

in which the terms can be seen to be closely analogous to the volume integral terms in 
(6). Thus using the expression (12) for the total momenta we can combine (6) and (55) 
so as to obtain (with the aid of (52) and (6.320) 1 ) the final result 

5M-0 H 5J BH -$ H 5Q H -^-5A = J k-i L x 5(n x -dT L )+ J (5tt x )- (kAn x ) • dE (56) 

(in which J BH = J H + J FH = J - J MC ). 

5. Interpretation in the generic (convective) case 

The first term on the right hand side of (56) can be interpreted as representing 
the static injection energy: it has the form 

J k-i L x 5(n x -dT L ) = J £ x 5(dN x ) (57) 

where diV x , as given by (42), is the flux of the xth constituent across dE, and where 
£ x is the canonical energy per particle, as given by (43) , which satisfies the conservation 
law 

n |x| -W*=0 (58) 

for each x- 

The form of the second term in (56) (which vanishes altogether in a static sit- 
uation, i.e. one with all current vectors parallel to k) may seem less familiar, but it 
is in fact a generalisation of the Faraday term IbT representing the energy needed to 
make a change b~T in the magnetic flux T through a circuit carrying a current I. The 
relationship can be made rather more apparent if we decompose each flux n x into a 
static part plus a spacelike part c x that is confined within the spacelike hypersurface 
E, in the form 

n x = nfk + c x , c x ■ dE = (59) 

since we may then write 

J (Stl x )- (kAn x ) -dE = -J (c x 5E. x )(k-dT L ) . (60) 
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We are always at liberty to decompose the bypersurface element dE in the form 

dE = dx- dS (61) 

where dx is an arbitrary infinitesimal displacement and dS_ is a normal 2 surface element 
transverse to the displacement. If we take dx to be locally parallel to one of the fluxes, 
c x say, we shall have 

-(c ]x] 6z x )(k-dE) = (dx-5n x )dC lxl (62) 

where dC |x| represents the current of the xth constituent through the element dS as 
defined by 

dC x = c x • (jfe- dS) . (63) 

Thus dC x represents the number of particles crossing the element dS_ per unit of the 
ignorable (not proper) group time t. In general one can envisage that the flux lines in E 
will have complicated ergodic properties. However if their global structure is such that 
they form sufficiently well behaved simple closed circuits then one can see from (62) 
that we shall have 



5M - H 5J BH - $ H 5Q H - ^SA = J £ x 5(dN x ) + J (5C x )dC^ 



(64) 



where for each circuit the conserved circulation integral C x , which is a generalisation of 
the Faraday flux, is as defined by (24). 

6. The circular (non - convective) case 

Despite the restriction imposed by the non-ergodicity hypothesis, the deceptively 
simple form (64) still represents quite a powerful generalisation compared not only with 
the original variation law of Bardeen et al (1973) 2 ' 3 , but even with the electromagnetic 
extension of Carter (1973b) 4 , since these earlier versions applied only to flows with 
strictly (not just topologically) circular circuits. To show how these earlier versions 
are obtained in the appropriate, non - convective (and non - conducting) limit, we first 
introduce a vector v x defined by 

c >»fVi ( 65 ) 

and representing the mean displacement of the particles of the xth constituent per unit 
of the ignorable (not proper) group time t. (If there were any particles in the immediate 
neighbourhood of the black hole, their 3-velocities as thus defined would evidently have 
to tend towards the Damour 7 velocity v H - as defined in section 6.3 1 - in the limit on 
the horizon.) We see that in terms of the corresponding unit flow vectors u x , defined 
by 

n x = n x u lx{ , |-u |x| | 2 = -1 (66) 

the 3-velocity 4-vectors will be given by 

k + v x =f {x{ u x (67) 
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where f x is the relevant time dilation (or "mean redshift") factor. We can use these 
3-velocities to regroup the terms on the right hand side of (56) so as to obtain the 
variation law 



5M 
where 



Q»SJ BH -^SQ H -^-SA = J** o) 6(dN x ) - Jj b x 6(7r x b n«dX a ) (68) 



n * :=-(k + v w )-Z x (69) 



which is valid quite generally (regardless of the topological behaviour of the flow lines) . 

If we now impose the restriction that the flow lines be strictly circular, so that 
the 3-velocities take the form 

v x =n x m (70) 
where Vt x are the corresponding angular velocities then we shall have 

5M - Sl*6J BH - $h 5 q h _ JLsA = £ n* 0) 6{dN x ) + £ O x 5(dJ M *) (71) 

with 

d^=^'dS (72) 

where the constituent angular momentum fluxes are given by (48) . The quantity 7r^ Q) , 
which is a generalisation of what Bardeen 2 has called the zero angular momentum 
injection energy, is given by 

' r ?o,=^- n m^ x =/*fo,+ eX *ix K o, (73) 

where the effective chemical potential /j,* is related to the corresponding local rest 
frame chemical potential (i x by 

=-(£ + tf |x| ).£* =f|x|/4 x (74) 

with 

V x = -U\x\-f± x ( 75 ) 
and where the corotating electric potential is given by 

$*(o) = -(£+v,M = ( 76 ) 

with 

^ x = -u lxl -A. (77) 
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